MATRICES RELATED TO DIRICHLET SERIES 
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DAVID A. CARDON 

Abstract. We attach a certain n x n matrix A„ to the Dirichlet series 
L{s) = X^fcLi o.kk^'^ ■ We study the determinant, characteristic polynomial, 
eigenvalues, and eigenvectors of these matrices. The determinant of An can 
be understood as a weighted sum of the first n coefficients of the Dirichlet 
series L{s)^^ . We give an interpretation of the partial sum of a Dirichlet 
series as a product of eigenvalues. In a special case, the determinant of A„ 
is the sum of the Mobius function. We disprove a conjecture of Barrett and 
Jarvis regarding the eigenvalues of ■ 

1. Introduction 
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To the Dirichlet series 



we attach the n x n matrix 



Lis) = E 



k=l 



k=l 

where En{k) is the nxn matrix whose ijth entry is 1 if j = ki and otherwise. 
For example, 

/oi 02 as 04 05 aeX 
fli 02 as 
ai a2 
ai 

ai 

V aj 

Since 

E^{k^)Er,{k2) = E„{k,k2) 

for every ki,k2 G N, formally manipulating linear combinations of En{k) is 
very similar to formally manipulating Dirichlet series. However, because En{k) 
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is the zero matrix whenever k > n, the sum defining D„ is guaranteed to 
converge. Of course, the n x n matrix contains less information than the 
Dirichlet series. Letting n tend to infinity produces semi-infinite matrices, the 
formal manipulation of which is exactly equivalent to formally manipulating 
Dirichlet series. 

Let Wn be the matrix whose first column is the weight vector (0, W2, W3, . . . , 
and whose other entries are zeros. Define the nxn matrix An (and the special 
cases Bn and C„) by 



(1) 



An = Wn + Dn, 

Bn = Wn + Dn whcu Ofc = 1 for all k, 

Cn = Wn + Dn whcu ttfc = 1 and Wk = 1 for all k. 



For example, Aq, Bq, and Cq are the following three matrices: 



/ ai a2 03 04 05 agX 

W2 ai a2 as 
W3 ai a2 
W4 ai 

\w6 ai/ 



(\ 1 1 1 1 IN^ 

W2 1 1 1 

W3 1 1 

W4 1 
W5 1 
V^6 1/ 



(\ \ \ \ \ \\ 
11 1 1 

1 1 1 
1 1 
1 1 

M 1/ 



We will always assume that a\ = \ since this ensures that the Dirichlet series 
^ akk'"^ has a formal inverse and since this is true for many Dirichlet series 
that arise in number theory. For notational convenience, we set Wi = 1, and 
occasionally we will write a{i) instead of Oj. Several authors have studied the 
matrices Bn and In [4j, it was observed that that 



(2) 



deti?n = Wfc/i(fc), 



k=l 



where fi is the Mobius /i-function. This is a special case of the slightly more 
general fact (see Theorem 12.11 below) that 



(3) 



det An = ^ Wkbk 



k=l 
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where the numbers are the coefficients of the formal series 

k=l 

Thus, det A„ is a weighted sum of the coefficients of L(s)~^. 

To obtain ([2]) from ([3]), choose the Dirichlet series to be the Riemann zeta 
function C^{s) = YlT=i ^'^ that = 1 for all k. This corresponds to the 
case of the matrix i?„. Since C{s)~^ = Xlfcli /^(^)^~'^' where /i is the Mobius 
/i-function, it follows that 6^ = fi{k). One particularly intriguing choice for Wk 
is Wk = k~''. Then ([3]) results in the truncated Dirichlet series 

k=l 

As the asymptotic growth of sums of the type in equation ([3]) is important 
to analytic number theory, representing those sums in terms of determinants 
becomes very interesting. 

Recall that the Riemann hypothesis is equivalent to the statement 

n 
k=l 

for every positive e. Thus, the Riemann hypothesis is equivalent to 

detC„ = 0{n^/^+'), 

for every positive e. 

In [1], Barrett, Forcade, and PoUington expressed the characteristic polyno- 
mial of Cn as 

where r = [\og2 n\ and where the coefficients v{n, k) were described in terms 
of directed graphs. We will refer to the eigenvalue 1, whose multiplicity is 
n — r — 1, as the trivial eigenvalue. The eigenvalues A 7^ 1 will be called 
nontrivial eigenvalues. In Theorem 13.21 we extend this result by determining 
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the characteristic polynomial of the more general matrix An- In jT], it was 
shown that the spectral radius p{Cn) of Cn is asymptotic to ^Jn. 

Barrett and Robinson [5] determined that the sizes of the Jordan blocks of 
Bn corresponding to the trivial eigenvalue 1 are 

Llog2(n/3)J + 1, Llog2(n/5)J + 1, . . . , Llog2(n/{n})J + 1, 

where {n} denotes the greatest odd integer < n. Theorem 14.11 of this paper 
shows that each nontrivial eigenvalue of An is simple and expresses a basis for 
the one-dimensional eigenspace in terms of a recursion involving the coefficients 
of Prn{x) for m < n, enhancing our understanding of the Jordan form of An- 
Theorem 14.21 gives a similar result for the transpose of An- 

The coefficients of the characteristic polynomial of Cn are related to the 
Riemann zeta function as follows: If {({s) — 1)*^ is expressed as a Dirichlet 

v-~^cx) c(m,k) J 1 1 

series > . ^ ' ' so that 

oo oo / oo / , X ' 

V^'V / / k=0 k=0 \m=l 



then 



v{n,k) = ^c{j,k). 



Evaluating Pn{x) at x = gives the fundamental relationship 

n Llog2 "J 

detC„ = 5^/i(z)= n ^= £ i-l)^vin,k), 

i=l A nontrivial k=Q 

where v{n, 0) is defined to equal 1. 

Barrett and Jarvis [2] showed that C„ has two large real eigenvalues A-t 
satisfying 

and that the remaining [log2 n\ — 1 small nontrivial eigenvalues satisfy 

|A| < log2_,n 
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for any small positive e and sufficiently large n. Based on numerical evidence 
for various values of n as large as n = 10^, they also made the following 
two-part conjecture: 

Conjecture 1.1 (Barrett and Jarvis [2]). The small nontrivial eigenvalues A 
ofCn satisfy 

(i) |A| < 1, and 
ill) Re(A) < 1. 

The statement Re(A) < 1 is, of course, weaker than the statement |A| < 1. 
Vaughan [6] refined the asymptotic formula ([5]) for the two large eigenvalues 
and showed, unconditionally, that the small eigenvalues satisfy 

|A| < {\ognf'\ 

and, upon the Riemann hypothesis, that the small eigenvalues satisfy 
(6) |A| <loglog(2 + n). 

He later showed [7] that C„ has nontrivial eigenvalues arbitrarily close to 1 for 
sufficiently large n, suggesting that a proof of Conjecture 11.11 would likely be 
quite subtle. 

Investigations of the Redheffer matrix have been extended to group theory 
by Humphries [3J and to partially ordered sets by Wilf [8]. 

In Theorem 15.31 we resolve Conjecture 11.11 by showing that both parts are 
false. There exist values of n for which a small eigenvalue A satisfies both 
|A| > 1 and Re(A) > 1. To accomplish this we computed the characteristic 
polynomials for An for values of n as large as n = 2'^^, which we describe in [|5l 

2. The determinant of An 
We now find the determinant of An- 

Theorem 2.1. Let Dn be the Dirichlet matrix associated with the formal 
Dirichlet series L{s) = Yl^=i^kk~'^ where ai = 1, and write L{s)~^ = YlT=i^kk 
Let Wn be the matrix whose first column is (0,^2, . . . ,Wn)^ and whose other 
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entries are zero. Let An = Wn + Dn as in ([T]). Also, let An = Wn + . 
Then 

n n 

(7) det{An) = ^Wkbk and det(y4„) = WfcQfc- 

k=l k=l 

Corollary 2.2. The choice Wk = 1 produces partial sums of coefficients of 
Dirichlet series: 

n n 

(8) det An = '^^bk and detA„ = ^^afc. 

k=l k=l 

Corollary 2.3. The choice Wk = gives truncations of the Dirichlet series 
L{s)~^ and L{s): 

^ b " 

(9) detA. = ^^ and deti. = ^^. 

fc=l k=l 

If s is a complex number at which L(s) and L{s)~^ converge, 
lim det A„ = and lim det A„ = L(s). 

n^oo n— >oo 

So, Corollary 12.31 says that we may interpret det An and det An as approxi- 
mating values of Dirichlet series. Since the determinant is the product of the 
eigenvalues, this relates values of Dirichlet series with eigenvalues of matrices. 

Proof of Theorem \2.1[ This is essentially the same argument as the one given 
in Redheffer's note [4J where he found the determinant of Bn- Since Dn is 
upper triangular with diagonal entry 1, det D„ = detD~^ = 1. Then 

det An = det D'^ det A„ = det D'^ det(Vr„ + Dn) = det{D-^Wn + /„). 

The matrix D~^Wn has zeros in columns 2 through n, and its (1, l)-entry 
is Ylk=2'^kbk- Thus, det An equals the (1,1) entry of D'^Wn + In which is 
Ylk=i '^kbk- Replacing Dn with in the argument gives det An = Ylk=i '^kO'k- 

□ 



matrices related to dirichlet series 7 

3. The characteristic polynomial of A„ 

The characteristic polynomial = det{InX — An) plays a significant role. 
Previously, p„(x) was obtained for the special case Cn in and |6j. In this 
section, we will determine the characteristic polynomial of the more general 
matrix An. The following definition will be instrumental in describing both 
the characteristic polynomial of An and its eigenvectors. 

Definition. For integers n > 1 and > 0, we define d{n, k) to be the Dirichlet 
series coefficients of {L{s) — 1)^ . That is, 

(10) (iw-i)'-ff:^)'-E%^- 

\fc=2 / n=l 

We define f (n, k) and v^in, k) to be the weighted sums: 

v{n, k) = ^w{j)d{j, k), and 

(11) ^ 

vi{n, k) = 2_^w{ji)d{j, k). 

Several cases of this definition are important to keep in mind: d{l,0) = 1 
and d{n,0) = for n > 1; also, both d{n,k) and v{n,k) are zero if n < 2*^ 
since a number smaller than 2'^ cannot be written as a product of k nontrivial 
factors. 

From the definition of d{n, k), 

d{n, k) _ fy^ f ^\ - ^\ ( c^l^, k - 1] 

n=l \fc=2 / \fc=2 / \fc=2 / \n=l 

which immediately gives the elementary recurrence relation: 
Lemma 3.1. Ifk>l, then 



[12) d{n, k) = a{i)d{n/i, k — 1) = a{n/j)d{j, k — 1). 



i\n j\n 
Ki j<n 



Theorem 3.2. The characteristic polynomial Pn{x) = det(x/„ — An) is 

r 

(13) Vn{x) = {x- \Y-'-^{{x - 1)'+^ -^v{n,U){x - \)''-^^^ 



k=l 
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where r — [log2nJ. Consequently, ifv{n,r) ^ 0, the algebraic multiplicity of 
the trivial eigenvalue X — 1 is n — r — 1. 

Proof. We will use the cofactor expansion to calculate the characteristic poly- 
nomial Pn{x) — det{xln — An). Write M„ — xin — and let 

Mn{ii,...,is I ji,...,jt) 

denoted the matrix obtained by removing the rows indexed by ii, . . . ,is and 
the columns indexed by ji, . . . ,jt from M„. The cofactor expansion of the 
determinant along the first column is 

n 

Pn{x) = (a;-l)" + 5^(-l)^WfcdetM,(A; | 1). 

k=2 

The matrix Mn{k | 1) is a block matrix whose upper left (A; — 1) x (A; — 1) 
block is Mk{k \ 1), whose lower left {n — k) x {k — 1) block is zero, and whose 
lower right {n — k) x {n — k) block is upper triangular with diagonal entries 
x-l. Thus 

detM„(A; | 1) = {x - lY'^ del Mk{k \ 1), 
where we understand detMi(l | 1) to be 1, and 

n 

(14) pn{x) = {x- If + Y,i-lfM^ - 1)""' det Mk{k I 1). 

k=2 

The £th entry in the last column of Mk{k \ 1) is —an/e if i divides k; oth- 
erwise, it is zero. Then the cofactor expansion of det Mk{k \ 1) along the last 
column is 

det Mk(k I 1) = ^(-l)'=+^afe/,detMfc(^, A; I l,k). 

e\k 

i<k 

The matrix Mj.{l, k\l,k) \s also a block matrix. Since the upper left (£ — 1) x 
- 1) block is Mi{l I 1), the lower left (fc - £ - 1) x (£ - 1) block is zero, and 
the lower right (A; — £ — 1) x {k — l—l) block is upper triangular with diagonal 
entries x — 

del Mk{l,k\ l,k) = {x - l^'^-^ det Me{i \ 1). 
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This shows that 



(15) det Mkik I 1) = - 1)'^-^-! det M^(£ | 1). 

e\k 

In other words, the quantity qk{x) = {—l)''^^Mk{k \ 1) satisfies the recurrence 
relation: 



qi{x) = 1, 

^^^^ qkix) = J2<^k/i{x-l)'-'-%{x) foTk>l. 

e\k 

i<k 

On the other hand, consider the polynomial defined by 



(17) U{x) = J2d{i,j){x-lY-^-\ 

Then ti(x) = 1. For i > 1, the term in the sum corresponding to j = is 
zero since d{i, 0) = in that case. For k > 1, calculating the right hand side 
of ffTBl) with ti{x) in place of qe{x) and applying Lemma [XT] gives 



Y,(^k/iix-l)'-'-\{x) = 5^5^a,/,rf(£,j)(x- 1)'="^-' 

e\k e\k j>o 

e<k e<k 

= J2d{k, J + l)ix-ir-^-' 

= Y,d{k,3){x-lf-'-' 
i>i 

= tk{x). 



Since t^i^x) and qk{x) both satisfy the same recurrence relations, they are equal. 
This shows that 



{-lf-^Mf,{k I 1) = g,(x) = tfc(x) = d{k,j){x - l 

i>o 

Substituting the last expression into (fT4|) gives 

n 

p„(x) = (x - 1)" - ^ 5^^,rf(A;, j)(x - ir-^-' 
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= (x-ir- J]t;(n,j)(a:-ir-^-^ 

i>i 

Since v{n,j) = for j > r = [log2(n)J, this is 

r 

r 

= {x- 1)"--^ ((x - 1)^+1 - J] ^(n, j) (x - l)'--^) , 

i=i 

which proves the theorem. □ 

4. The eigenvectors of An 

Theorem 4.1. Let X 1 be a nontrivial eigenvalue of An. Then X is a simple 
eigenvalue, and a basis for the one dimensional eigenspace of An associated 
with X is the vector 

u= [X-l,X2{[n/2\),Xs{[n/3\),X,{[n/4\),...,Xn{[n/n\)f 

where 



Proof of Theorem 4-T For i > 2, the ith entry of y4„M is 



{AnU)i = Wi{X - 1) + ^ a£M« 

l<£<n/j 

= w;,(A-l)+ 5^ aeXu{[n/{a)\) 

l<l<n/i 



W 



,(A-1)+ J2 Yl wii^'rn)d{m,k){X-l)-^ 

l<i<n/i k>0 

l<m<n/i 

(A-l) + ^( J2 a{i)w{iim)d{m,k)yX-l)-'' 

k>0 l<i<n/i 
l<m<n/{ei) 

(A-l) + ^( J2 w{tt)J2ait/s)d{s,k)yX-iy'' [sett 



k>0 ' l<t<n/i s\t 

WiiX-l) + J2 Yl ^^i^t){dit,k) + d{t,k + l)){X-l)-'' [by (HI 

k>0 l<t<n/i 
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fc>0 fc>l 
fc>0 fc>0 

= XX,{[n/t\) 
= Xui. 

In the calculation for (A„-u)j with i > 2, the term ajW^j when i = 1 was equal to 
aiXi{\n/i\), but this term should be omitted from the case i = 1. Taking this 
into account and going to the second to last step of the previous calculation 
gives 

= (^-i)(i + E^(^'^)(^-i)"') 

k>l 

= (A - 1) [1 + (A - 1)] by Theorem [3J 
= A(A-1) 
= Ami- 

This shows that the vector m is a nonzero eigenvector for A. 

To see why the eigenspace of A is one-dimensional, consider the submatrix 
of An — XI obtained by deleting the first row and column. This (n — 1) x (n — 1) 
matrix is upper triangular with nonzero entries on the diagonal. Hence, it is 
invertible implying that the rank of An — XIn is > n — 1. Since we found a 
nontrivial eigenvector, the nullity is > 1. So, the nullity of An — XI must be 
exactly one. This completes the proof. □ 

Theorem 4.2. Let X ^ 1 be a nontrivial eigenvalue of An- A basis for the 
one dimensional eigenspace of A^ associated with X is the vector 

(18) v=[l,Y,{2),Y,{3),...,Y,{n)f. 
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where 



^ d{q,k) d{q,l) d{q,2) 

Interestingly, the algebraic expression for v does not explicitly rely on the 
symbols W2, . . . ,Wn in the first column of An- However, altering W2, . . . ,Wn 
changes the possible numeric values of A. 

Proof of Theorem \4.^ For i > 2,the ith entry of A^v is 

i\i 

fc>0 i\i 

= ^[d(^,fc) + rf(^,fc + l)](A-l)-^ by ([12]) 

fc>0 



F(z) + (A-l)5^d(z,A;)(A-i; 



k>\ 



Y{i) + (A - \)Y{%) [since d{i, 0) = 0] 
\Y{i). 



The first entry of Af^v is 



l<i<n 
k>0 l<j<n 



k>0 



k>l 



1 + (A - 1) [by Theorem [3J] 
Xvi. 
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This shows that v = [Y{1), . . . ,Y{n)]'^ is a nonzero eigenvector of A^. The 
dimension of the eigenspace is one, as explained in the proof of Theorem 14.11 

□ 

5. Computing eigenvalues of C„ for large n 

Theorem 13.21 expresses the characteristic polynomial of the matrix An in 
terms of the numbers v{n, k). In this section, we will explain how to explicitly 
calculate the characteristic polynomial Pn{x) for large values of n for the special 
case Cn in which = = 1 for all i. The method given below in Theorem l5.2l 
was used to find Pn{x) for n as large as n = 2^^ in a few hours on a desktop 
computer. To accomplish this, it is necessary to use a more efficient algorithm 
for finding the coefficients than a brute force approach based directly on the 
definition of matrix C„. Even with Theorem 13.21 we need a better method 
for computing v{n,k) than the direct application of the definition of v{n,k) 
in (mi). 

Lemma 5.1. Suppose ae = W£ = 1 for all i. // 1 < 2^ < n, then 

(19) v{n, k) = J2v{[^\,k-l)=J2{[j\- IwilH^, k - 1). 

If both Ofc = 1 and Wk = 1 for all k, then v{n, k) represents the number 
of ways to form products of k nontrivial factors whose product is < n and 
where order matters. In this case, f (n, k) represents a count of lattices points 
in fc-dimensional space: 

(20) v{n,k) = |{(£i,...,4) : ■ ■ ■ 4 < « and £i > 2 for alH}|- 



Proof. The first equality in ( 1191) is evident from (12011 by letting one component 
of (£i, . . . , £fc), say be the index of summation i. The second equahty in (fT9l) 
is obtained by re-indexing the sum over the distinct values of j = \n/i\. For 
a given positive integer j, 



n n n 

j<-<j + l -^^ <i < -. 

« J + 1 J 



3 = 

Thus, the number of distinct i for which \n/i\ = j is ["J — L^^J- ^ 
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The first recursion formula in f|T9l) is computationally inefficient since there 
can be many distinct values of ii and 22 such that = [n/i2\ ■ The second 

is inefficient since there can be many values of j such that [n/j\ — [n/{j + 1)J 
is zero. The next theorem helps to remove this redundancy by rewriting the 
summation to have significantly fewer terms. 

Theorem 5.2. Assume ai = wi = 1 for all i. Suppose 1 < 2^ < n and k > 1. 

Then 

Lv^J 

(21) vin,k) = J2<m,k-l)+ E {[j\-[jTT\)^(3'k-l), 

i=2 j=2fe-i 

Proof. This argument applies the hyperbola method from analytic number 
theory. Rewrite ffTIJl) as 

(22) v{n,k)= J2 ^([fJ'^-l)+ E 

Ln/jJ>Lv^J+l ln/i\<lV^\ 

where the index i in each summation satisfies 2 < i < [n/2'^^^\. In the first 
summation, since both i and [y/n\ + 1 are integers, 



n 



I 



n . ^. n 



> [v^J + 1 ^ - > Lv^J +1 ^ i< ^ ^ ^ i < 



n 



This gives the value s = \n/{ [v^J +1)J iii the first summation in equation (12T|) . 
In the second summation in fl22l) . we re-index the sum over the distinct values 
of j = \ n/i\ < [v^J- For a given positive integer j, 



J 



n n n 

^ j<-<j + l 4^ <i<-. 

I j + 1 j 



Thus, the number of distinct i for which \n/i\ = j is [yj — Lj^J' allowing 
the second summation in (l22l) to be written as 



E ([fJ-lTTTJ)-^^-!)- 



This proves (ETil . □ 
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It is interesting to note that s in Lemma 15.21 is equal to either [v^J or 
IV^\ ~ 1 according to 



n 



LLv^J + 1, 



Lv^J iin-[^P>[^\., 
[v^J - 1 if n - [v^j2 < [^J. 



Theorem 5.3. Conjecture \1.1\ is false. There exist values of n for which a 
small eigenvalue A o/C„ satisfies both |A| > 1 and Re(A) > 1. 



Proof. The characteristic polynomial Pn{x) of the general matrix An was given 
in Theorem I3.2[ By implementing the recursive formula in Theorem I5.2[ we 
were able to calculate the characteristic polynomial for the special case C„ for 
relatively large values of n, such as = 2^^, within a few hours on a desktop 
computer. 

A table listing the maximum absolute value and real part of small nontrivial 
eigenvalues of C„ for n = 10^ and n = 2^ with 28 < r < 36 is given below: 



n 


max{ A } 


max{Re(A)} 




1,000,000 


0.983108 


0.983108 


228 


268,435,456 


0.998885 


0.998739 


229 


536,870,912 


0.999120 


0.998989 


230 


= 1,073,741,824 


0.999324 


0.999206 


231 


= 2, 147, 483, 648 


0.999501 


0.999395 


232 


= 4, 294, 967, 296 


0.999676 


0.999560 


233 


= 8,589,934,592 


1.002646 


0.999704 


234 


= 17,179,869,184 


1.005213 


0.999829 


235 


= 34,359,738,368 


1.007423 


0.999939 


236 


= 68,719,476,736 


1.031192 


1.000036 



The example with n = 2^^ provides a counter-example to both parts of Con- 
jecture [LTI A sample of the coefficients v{n, k) of Pn{x) for n = 10®, n = 2^^, 
and n = 2^® is given in Table [H □ 
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Table 1. Values of v{n, k) for n = 10^ 



k 



n = 2^8, and n = 2 



i;(236,A;) 



36 



1 

2 

3 
4 
5 
6 
7 
8 
9 

10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 



999999 

11970035 
67120491 
233959922 
567345854 
1015020739 
1386286166 
1475169888 
1237295133 
822451796 
433656192 
180821164 
59146673 
14935574 
2829114 
383693 
34630 
1672 
20 



268435455 

4714411991 
39550266080 
210866000001 
801946179797 
2314766752399 
5267935378357 
9693670870002 
14675212443928 
18500845515388 
19585798031078 
17506983509953 
13254336924806 
8508754910066 
4628591443629 
2128656115076 
824357770148 
267263904116 
71941723387 
15889930335 
2830811858 
396537923 
42162106 
3284753 
177731 
4707 
55 
1 



68719476735 

1587951104025 
17712699735807 
127006997038631 
657738684402616 
2620541404211325 
8354699452581663 
21888970237054221 
48028484118248949 
89496511738284187 
143118705146069804 
197979547265239162 
238336089820847725 
250812663743567239 
231467885026020936 
187727209728498411 
133949812310943213 
84103735312636462 
46433832280215021 
22505741596654059 
9551600816612963 
3536981261202340 
1137490727898326 
315879734318303 
75228001661856 
15244074212812 
2604780031507 
371154513760 
43388420848 
4049932603 
290175811 
15487073 
582143 
9555 
71 
1 
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